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ABSTRACT 

Aims. 4U 1954+319 was discovered 25 years ago, but only recently has a clear picture of its nature begun to emerge. We present 
for the first time a broad-band spectrum of the source and a detailed timing study using more than one year of monitoring 
data. 

Methods. The timing and spectral analysis was done using publicly available Swift, INTEGRAL, BeppoSAX, and RXTE /ASM 
data in the 0.7 to 150 keV energy band. 

Results. The source spectrum is described well by a highly absorbed (Nh ~10 23 cm -2 ) power law with a high-energy 
exponential cutoff around 15 keV. An additional black body component is needed below 3 keV to account for a soft excess. 
The derived ~5 hr periodicity, with a spin-up timescale of ~ 25 years, could be identified as the neutron star spin period. The 
spectral and timing characteristics indicate that we are dealing both with the slowest established wind-accreting X-ray pulsar 
and with the second confirmed member of the emerging class dubbed "symbiotic low mass X-ray binaries" to host a neutron star. 
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1. Introduction 

Since its discovery by the ARIEL satellite in 
1981, 4U 1954 +319 has appeared as a flaring 
X-ray source (jWarwick et al. 198l[) . Located in 
the Cygnus region, this source was observed by 
different X-ray telescope s llForman et al. 1978 ; 
Cook. Warwick fc Watson 19841: IWarwick et al. 1988 : 



Tweedy et al. 19891 IVoees et al. 1999|) . From the ob 
served spectral behaviour, the source was tentatively 
classified as a high-mass X-ray binary system hosting 
a neutron star (NS). However, little was known about 
the system, as the optical counterpart was not identi- 
fied due to the ambiguous determination of the source 
position. The companion star was then classified as a 
close M-type giant or a distant and reddened Be star 
( Cook. Warwick fc Watson 1984HTweedv et al. 1989[) . 



Recently, u sing t he Chandra refined source position, 
iMasetti et all l|2006j) used optical spectroscopy to iden- 
tify its companion as an M4-5 III star, which is located 
within 1.7 kpc. Such identification suggests the system 
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is composed of a compact object accreting through the 
wind of its M-type giant companion star. This makes 
4U 1954+319 the third low mass X-ray binary (LMXB) 
hosting an NS and a late type giant co mpanion, after 
GX 1+4 (e.g . _ IChakrabartv fc Roche 19971) and possibly 
4U 1700+24 <|Gaudenzi fc Polcaro 1999|) . for which how- 
ever no coherent pulsat ion has been reported to date 
(e.g. IMasetti et al. 2002f) . Therefore 4U 1954+319 could 
be attributed to the so-called "symbiotic X-ray binaries" , 
an emerging subclass of LMXBs. These systems with an 
evolved giant donor are the probable progenitor s of most 
wide-orbit LMXBs l|Chakrabartv fc Roche 1997|) . 



Only recently has a periodic signal been reported 
from this source, whit a period of ~ 5.09 h ours quasi- 
monotonically decreasing l|Corbet et al. 200rih . If this pe- 
riod can be attributed to the spin, 4U 1954+319 would be 
the slowest binary NS known, with the possible exception 
of the enigmatic sou rce IE 161348-5055 i n the Supernova 
remnant RCW 103 l|De Luca et al. 20061) . In this Letter 
we used several high-energy observations to unveil the na- 
ture of 4U 1954+319 through spectroscopic and timing 
analysis. 



2 



F. Mattana, D. Gotz, M. Falanga, et al.: 4U 1954+319 



2. Observations and data analysis 

The INTEGRAL data set was obtained using the 
2003/2005 publicly available observations within < 
12° from the source direction. We analysed 613 
pointings from the IBIS /ISGRI coded mask telescope 
ljUbertini et al. 20"03l iLebrun et al. 2003^ at energies be- 
tween 18 and 200 keV, for a total exposure time of 1.4 
Ms. The data reduction was performed using the standard 
Offline Science Analysis (OSA) version 5.1. We extracted 
the 20-60 keV band light curve based on the single point- 
ings. During the INTEGRAL observation, 4U 1954+319 
was mainly in a quiescent phase (52717-53144 MJD), and 
in the last part of the observation, it entered a strong 
outburst phase (53312-53708 MJD) with an increase in 
count rate of a factor 40. During the outburst phase, 
4U 1954+319 is clearly detected in the mosaic image with 
a significance level of ~ 45c. During the quiescent phase, 
4U 1954+319 was not detected at a statistically signifi- 
cant level in the ISGRI data. Therefore we excluded these 
data from the timing and spectral analyses. 

We also included for the spectral analysis the publicly 
available data from the Narrow Field Instruments (NFIs) 
on board BeppoSAX. The source 4U 1954+319 was ob- 
served on May 4, 1998 for a net exposure of 19 ks in the 
LECS llParmar et al. 1997^ and for 46 ks in the MECS 
(Bo ella et al. 19971) . The average spectrum of the source 
was extracted from a region of 6' in both instruments. 
The spectra were then rebinned in order to have at least 
30 counts per channel and three channels per resolution el- 
ement. Note that BeppoSAX observations were too short 
to enable the timing analysis. 

We also analysed the Swift/BAT 

1 Barthelmv et al. 2005h survey data, collected be- 
tween December 2004 and April 2006. The analysis was 
performed using the HEASoft software (release v. 6. 0.5). 
For each single pointing (with typical integration time 
of 5 minutes), we extracted an image in the 15-50 keV 
energy range and evaluated the source count rate at the 
Chandra position. Thus, we obtained a barycentered light 
curve containing ~ 23000 data points. 

In Fig. ^ we show the Rossi X-ray Timing 
Explorer (RXTE) All-Sky Monitor (ASM) light curve of 
4U 1954+319, where we report the time interval in which 
the INTEGRAL, Swift, and BeppoSAX observations were 
performed. 



2.1. Spectral analysis 

For the spectral analysis, we used ISGRI (18-150 keV), 
LECS (0.7-4.5 keV), and MECS (1.7-10 keV) data. The 
non-imaging BeppoSAX spectrometer PDS data (15-300 
keV) could not be used due to the contamination of Cyg 
X-l. The ana lysis was done using XSPEC version 11.3 
l)Arnaud 1996h . All spectral uncertainties in the results 
are given at a 90% confidence level for single parameters. 

We first computed a hardness-intensity diagram us- 
ing the RXTE/ ASM data. The hardness is the ratio of 
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Fig. 1. The 4U 1954+319 RXTE/ ASM light curve in the 
1.5-12 keV energy band (data-averaged over a 7-day in- 
terval). We conv erted the ASM coun t rate to flux using 1 
Crab w 75 cts/s l|Levine et al. 1996|) . 



the count rates in the 1.5-5 keV to the ones in the 5- 
12 keV energy band, and the intensity is the 1.5-12 keV 
count rate. Each point corresponds to ~ 7 days of inte- 
gration time. The source shows a significant linear corre- 
lation (slope w 1) between flux and hardness. Therefore 
we searched for spectral changes in the ISGRI data, di- 
viding the observation into four different intensity levels. 
We fitted each subset with a simple power-law (pl) model, 
but it turned out to be statistically inadequate (x 2 /d.o.f 
= 43/12). The best fit was found by replacing the PL with 
a cutoff PL model thereby obtaining a x 2 /d.o.f = 10/11. 
We found the spectral parameters for the four data sets 
to be the same within the error bars, most likely due to 
the low statistics during the flare, so we did not observe 
the same hardness-intensity correlation as observed below 
12 keV with the ASM data. The best-fit values for the 
average ISGRI spectrum are found for a PL photon index, 
T, of 0.7 ± 0.1 and a high-energy cutoff E c ~ 14±§ keV. 

We then fitted the BeppoSAX spectrum using a sim- 
ple photoelectrically-absorbed PL model plus a black body 
(bb) model for the soft excess below 3 keV, as oft en ob- 
served in accreting pulsars (e.g. iHickox et al. 200*3) and a 
Gaussian emission line, resulting in a x 2 /d.o.f. = 132/71. 
The best fit was found for a heavily absorbed spectrum 
with Ah = 2.8 ± 0.2 x 10 23 cm -2 , a bb temperature of 
fcT=0.11±0.01 keV to model the soft X-ray excess, a pl 
index of T = 1.5±0.5, E c ~ 17±12 keV, and an iron emis- 
sion line at ~6.33±0.15 keV (with the line width fixed at 
0.5 keV). The quite high reduced x 2 indicates that the 
soft excess is not perfectly modeled by the BB component. 
The hydrogen column density, Ah, is found to be 30 times 



high er than the Galactic valu e reported in the radio maps 



higher than the (jalaci 
of iDickev fc LockmarJ 



The cutoff value is not well-constrained using the NFIs 
or ISGRI alone. This is due to the fact that it falls outside 
of the energy ranges of the individual instruments. In order 
to better constrain the broad-band (0.6-200 keV) spectral 
characteristics we jointly fitted the NFI and the ISGRI 
data. A multiplicative factor for each instrument was in- 
cluded in the fit to account for the uncertainty in the cross- 
calibration of the instruments, as well as variability across 
the non-simultaneous observations. The best-fit parame- 
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Fig. 2. Broad-band spectrum of 4U 1954+319, fitted with 
an absorbed bb Gaussian line, and cutoff PL. The data 
points correspond to the LECS/MECS (0.7-10 keV) and 
ISGRI (18-150 keV) spectra. The lower panel shows the 
residuals with respect to the best-fit model. 



ters are T=l. 1+0.1, £ c =17±2 keV, fcT=0.123±0.004 keV, 
and a similar value of Nh, with a x 2 /d.o.f. of 138/83. The 
broad-band spectrum, together with the best-fit model, 
are shown in Fig. [21 

2.2. Timing analysis 

For the timing analysis we used the BAT (15-50 keV) 
and ISGRI (18-50 keV) light curves after solar-system 
barycentric correction. 

We searched for coherent pulsations of the source in the 
BAT data by computing a power density spectrum (PDS) 
in the frequency range between 2x 10 -8 and 10 Hz from 
fast Fourier transforms. In the resulting PDS, a signal is 
evident at v = 5.45 x 10~ 5 Hz. The peak is broad with an 
FWHM of 6.7xl0~ 7 Hz, suggesting a period evolution. 
Therefore the light curves were grouped into 100 cycles 
per time interval, and the best period was determined us- 
ing an epoch-folding analysis. The distribution of the y 2 
values versus trial period were fitted as described in lLeahvl 
l(l987l) . The resulting 14 best-period values clearly show a 
spin-up trend, see Fig. |21 

Using the ISGRI light curves, we confirmed the pres- 
ence of a period at ~5.17 hr. We attempted to search for 
a period change in the ISGRI data; but due to the sparse 
coverage using the 100 cycle time intervals, we could only 
derive 2 significant measurements (see Fig. |2Jl , which are 
consistent with the previous findings. 

Fitting the period evolution with a linear function, the 
derived spin- up trend is around — 3xl0~ 5 s s _1 , which 
implies a spin-up time scale (P/P) of ~ 25 years. However, 
we note that some scattering around this trend is present, 
indicating a possible local spin-down. On the other hand, 
a much better fit was obtained using a sinusoidal function, 
indicating a possible orbital period of 385±4 days and a 
projected semimajor axis, a x sini, of 178±4 AU. However, 
since this periodicity is of the order of the time interval 
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Fig. 3. Upper panel: BAT light curve (15-50 keV) of 
4U 1954+319. The time bin size is 22 days. Lower panel: 
period evolution as measured with BAT (crosses) and 
ISGRI (diamonds). The inset represents the normalized 
BAT phase folded light curve between 53400 and 53500 
MJD considering P and P. Two cycles are shown for clar- 
ity. 

spanned by the data, it is difficult to identify it with a 
possible orbital period (see below). 

3. Discussion 

The timing and broad-band spectrum of 4U 1954+319 al- 
lowed us to investigate for the first time the nature of 
this system, which only recently has been tentatively at- 
tributed to the emerging NS "symbiotic LMXBs" class. 

The ~5 hr measured periodicity could be due either 
to the system orbital modulation or to the rotation of the 
compact object. We can rule out that we are measuring 
the system orbital period. In fact, given the typical lumi - 
nosity and temperature of an M 4-5 III star l|Lang 1 992). 
the stellar e volutionary models with solar metallicities 
l)Glaret 2004h predict a mass of ~1.2 M Q for the compan- 
ion. Assuming a mass within 3 M© for the compact ob- 
ject, the orbital separation would be much smaller than 
the donor star radius (~80 R©). We are then left with 
the hypothesis of measuring the spin period of the NS 
hosted in the system (a white dwarf would be excluded 
due to the excessiv e torque required, as already noted by 
lOorbet et al. 200rfl . 

Furthermore, the observed spin-up could be attributed 
to the accretion torque on the magnetized NS. While sug- 
gestive of a sinusoidal trend (see Fig. and Sect. 12.2(1 . the 
time evolution of the period cannot be accounted for by 
the Doppler shift of the system's orbital motion: in this 
case the derived orbital parameters would imply a com- 
panion star mass of a few 10 6 M©. Thus we can consider ~ 
400 days a lower limit to the orbital period of the system. 
This implies a wide orbit (a >2x 10 13 cm) a nd explains the 
fact th at the optical spectrum reported bv lMasetti et alj 
(2006) shows no sign of the influence of the X-ray source. 

We found a bro ad-band spectrum typica l of accreting 
X-ray pulsars (e.g. Ijoss &, Rappaport 1984j) and, in par- 
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ticular, similar to the ones of GX 1+4 llPaul et al. 20051) . 
and 4U 1700+24 l|Tiengo et al. 2005|) . The high column 
density found of the order of 3xl0 23 cm -2 could be at- 
tributed to the local (circumstellar) absorption. In fact, 
this value is typical of the local environment of symbi- 
otic stars as indicated by observations o f Rayleigh scat- 
tering in their UV spectra IjSchmid 19971) . The measured 
luminosity of ~ 2 x 10 35 erg s , assuming a 1.7 kpc 
source distance, points towards a wind-fed accretion sys- 
tem with the companion not filling its Roche lobe. Taking 
the period, its derivative, and the luminosity into ac- 
cou nt, one can u se the standard accretion torque mod- 
els IjHenrics 1 983') to derive an order-of-magnitude value 
for the dipolar magnetic field of the NS. It turns out to 
be B r~j 10 12 G, which implies a magnetospheric radius 
r m = 3 x 10 8 cm (L x /10 37 erg s-^-^B/lO 12 G) 4 / 7 ~ 
10 9 cm, smaller than the corotation radius (r co = 
(GM x P 2 pm /47r 2 ) 1 / 3 - 1.2 x 10 11 cm, where M x ~1.4 M© 
is the mass of the NS), allowing the accretion process. 
We are hence dealing with a regular wind accreting X-ray 
pulsar but with the longest known spin period. 

The very long period measured in 4U 1954+319 can 
be explained by taking into account the age of the com- 
panion star (~9xl0 9 vr. lClaret 2Qol and the standard bi- 
nary pulsar evolutionary picture l|Davies fc Pringle 198lj) . 
In fact, a newborn NS experiences a spin-down initially 
due to the magnetic dipole braking (radio phase). In our 
case this phase lasts ~8xl0 9 yr, namely the time the 
donor spends on the main sequence. Since the orbit is 
quite wide, the NS can spin down practically as if it was 
isolated. When the donor star leaves the main sequence, 
starting to evolve rapidly into a red giant, the star mass 
loss increases and there is coupling between the infalling 
matter stopped by the mag netosphere and the magneto - 
sphere itself (propeller phase. liTiarionov fc Sunvaev 1975|) . 
By means of these processes during its whole lifetime, an 
NS with a magnetic field of a few 10 12 G can reach periods 
of the order of 10 4 s in -8.5xl0 9 yr. 

The distance where the ram pressure of the accretion 
flow balances the magnetic pressure is called the Alfven 
radius. A further increase in the accretion rate, as the 
donor becomes a red giant, allows the corotation radius 
to overcome the Alfven radius and the NS to start ac- 
creting, becoming observable in X-rays and spinning-up. 
The emerging picture is an old, hence slow, NS accreting 
from the slow and dense wind of an evolved M type gi- 
ant on a wide orbit. The inhomogeneities of the wind can 
explain the observed long-term X-ray variability (see Fig. 
nj . The wide orbit also allowed the companion star not to 
be affected by the Supernova explosion that generated the 
NS thereby enabling the donor star to follow its natural 
evolutionary track. 

4. Conclusions 

We analysed most of the recently available high-energy 
data of 4U 1954+319 deriving for the first time its broad- 
band (0.7-150 keV) spectrum and analysing its timing 



characteristics. The luminosity and spin period derivative 
of the system allowed us to infer a pulsar magnetic field of 
the order of 10 12 G. In the framework of standard binary 
pulsar evolutionary models, this value, coupled to the age 
of the companion star, can justify the very long observed 
NS period without invoking any additional spin-down pro- 
cess. 

It turns out that 4U 1954+319 is the slowest accretion 
powered NS, and it is hosted in a "symbiotic LMXB" sys- 
tem. These are rare systems that could eventually evolve 
into wide-orbit binary systems hosting an NS and a white 
dwarf, but without mass transfer so hence very difficult to 
observe. 
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